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Neutrino telescopes with cubic kilometer volume have the potential to discover new particles. 
Among them are next to lightest supersymmetric (NLSPs) and next to lightest Kaluza-Klein 
(NLKPs) particles. Two NLSPs or NLKPs will transverse the detector simultaneously produc- 
ing parallel charged tracks. The track separation inside the detector can be a few hundred meters. 
As these particles might propagate a few thousand kilometers before reaching the detector, multiple 
scattering could enhance the pair separation at the detector. We find that the multiple scattering 
will alter the separation distribution enough to increase the number of NLKP pairs separated by 
more than 100 meters (a reasonable experimental cut) by up to 46% depending on the NLKP mass. 
Vertical upcoming NLSPs will have their separation increased by 24% due to multiple scattering. 



PACS numbers: ll.30.pb, 13.15-|-g, 12.60.jv, 95.30.Cq 

I. INTRODUCTION 

One of the contemporary questions in particle physics 
regards extensions of the standard model at higher en- 
ergy scales. New physics should appear around TeV en- 
ergy scales. Large telescopes for high-energy astrophys- 
ical neutrinos can complement colliders in searching for 
signatures of extension models. In some models, their 
sensitivity can extend beyond the reach of the LHC. More 
specifically it was shown that neutrino telescopes can di- 
rectly detect next to lightest supersymmetric particles 
(NLSP) [H E] and next to lightest Kaluza-Klein (KK) 
particles (NLKP) [3]. In scenarios were the lightest su- 
persymmetric (KK) particle is the gravitino (KK mode of 
the graviton) , charged right-handed sleptons (KK lepton 
modes) will be the NLSP (NLKP). In the supersymmet- 
ric scenario, neutrino telescopes searches will complement 
collider searches in the determination of the supersym- 
metry breaking scale [T . 

Searches for both NLSPs and NLKPs complement dark 
matter searches. The lightest particle in each of these 
scenarios arc dark matter candidates. Both the grav- 
itino and the KK mode of the graviton cannot be de- 
tected, at least with the current technologies. Therefore 
NLSP (NLKP) searches also probe dark matter, albeit 
indirectly. Other supersymmetry scenarios with charged 
leptons as NLSPs to be probed by neutrino telescopes 
were considered by [4j . These assume scenarios where ei- 
ther the neutralino or a super weakly massive interacting 
particle is the LSP. 

High energy (above 10'^ GeV) neutrino interactions in 
the Earth may produce supersymmetric or KK particles. 
These will go through a decay chain [3j which will 
always lead to the production of two NLSPs or NLKPs. 
These next to lightest particles will propagate through 
the Earth, eventually reaching a neutrino telescope. Due 



to their heavy mass, their specific energy loss (dE/dx) as 
they propagate through the earth is much smaller \T, 'S] 
than the loss by standard leptons. For this reason NLSPs 
and NLKPs may have a range of thousands of kilometers, 
while even the most energetic muons have a range of a few 
tens of kilometers. The smaller energy loss compensates 
the small NLSP and NLKP production cross sections. 

One signature for these new particles is two charged 
tracks transversing the detector [51 [3]. Here we inves- 
tigate the NLSP and NLKP multiple scattering while 
propagating through the Earth. We consider the sim- 
ulated distribution of their production point and energy 
loss as described in |5] The separation distance be- 
tween the two charged tracks in the detector will be a 
powerful discriminating observable. It is an excellent pa- 
rameter to separate NSLPs or NLKPs from the main 
background, which consists of dimuon events. As a pre- 
view to our conclusion, we note that multiple scattering 
will increase the number of NKLP pairs separated by 
more than 100 meters up to 46% depending on its mass. 
Although the NLSP pair separation increase is smaller 
(3%) the multiple scatter increases the number of verti- 
cal upcoming NLSPs by 24%. 

In the next section we describe how we account for 
multiple scattering and energy loss. Then we show our 
results, followed by a discussion on the detectability of 
these particle. Our conclusions follow in the last section. 



II. MULTIPLE SCATTERING 

The formulae for multiple scattering are well known; 
the scattering is independent of the particle mass. The 
angular deviation by a singly charged relativistic particle 
with momentum p passing through a slab with thick- 
ness L made out of a material with radiation length Xq 
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mass < A >, density < p > and radiation length for each 
of the Earth's layers. The layers are as in the density profile 
given in [S] . The height corresponds to the distance from the 
center of the Earth to the beginning of the layer. Layer 1 
starts at the center of the Earth. The top of the lO*'' layer 
is at the surface of the Earth (6378 Km). The values for the 
crust, mantle and core are given respectively in [9l ll0lfTT| . 

in space (i.e., the total angular deviation, not that in a 
plane) is well described by a Gaussian distribution [6]. 
The scattering angle in a given plane 9q is defined as 

The radiation length, Xq is calculated using Dahl's fit 
given in Equation 27.22 of [6J to the data: 

^ ^ 716.4 g cm-^A 
° Z{Z +l)ln{287/VZ) 

where A and Z are the atomic mass and number of the 
target. This equation agrees with Tsai's [7] values for 
Earth targets to better than 2.5%. 

Since the earth is inhomogeneous, we use the density 
profile described in [8]. This profile divides the Earth in 
ten layers. The Earth structure, mean atomic number 
{Z) , mass (A) of each layer depend on if it belongs to the 
crust, the mantle or the core of the Earth. These together 
with the mean density (p) and the radiation length Xq 
are given in Table [T] 

We calculate NLSP/NLKP energy loss following [21[3], 
and assuming that it is smooth. The stochastic nature of 
the energy loss might broaden the distributions slightly, 
but should not be a significant factor here. 

We calculate the NLSP/NLKP divergence from the 
original neutrino trajectory. The deviation due to the 
opening angle of the particle pair has been previously 
considered. Here, we consider the additional deviation 
due to the multiple scattering of the NLSP pairs. 

For each layer, the mean scattering angle is calculated 
from Eq. |2.1[ and, from that angle, the deviation from 
the neutrino trajectory is projected to a surface detector. 
The average deviation from the neutrino trajectory by 



a NLSP/NKLP particle due to scattering in layer i is 
di = where D is the distance from the middle of 

the layer to the surface; the -\/2 accounts for scattering 
in the x and y directions perpendicular to the direction 
of travel. 

As in Ref. |2, we integrate over different production 
points in the earth. 

The deviation in each layer is treated independently. 
Since the angles are independent, the deviations are 
added in quadrature. As we neglect the small initial en- 
ergy difference between the two particles, and as they 
scatter independently, their deviations are also added in 
quadrature and determine the total separation due to 
multiple scattering. This separation is then added in 
quadrature (because the angles are uncorrelated) with 
the separation due to the initial opening angle. 

Three factors affect the contribution to the multiple 
scattering from each layer. The inner layers are the dens- 
est, and also have the longest lever arm, increasing their 
contribution. However, the NLSP also have the higher 
energies while they are near the production point, de- 
creasing the scattering. Overall, the inner layers still 
make the largest contribution to the scattering - the lever 
arm wins out. 

We also considered the effect of magnetic fields on the 
separation; since the two NLSP have opposite charges, 
they are bent in opposite directions by the earth's mag- 
netic field. However, even after accounting for the in- 
creased field strength deep in the earth [TT, this does 
not contribute significantly to the scattering. 



III. RESULTS 

The analysis of NLKP detection by neutrino obser- 
vatories [3j considered three mass values of 300, 600 
and 900 GeV. The two NLKP track separation in the 
detector is the only discriminating parameter for the 
highest mass considered although not the only one for 
NLKP masses below 600 GeV. The background consists 
of dimuon events, mostly from charm decay. 

The NLKP energy loss depends on its mass. The larger 
the mass the less its energy will degrade while the particle 
propagates through the Earth. This means that heavier 
NLKPs can be produced farther from the detector. On 
one hand this enhances the multiple scattering, since the 
particle will have a longer lever arm. However on the 
other hand the NLKP enters each layer with larger en- 
ergy, reducing the scattering. 

Figures [T] and |2] show the effect of multiple scattering 
on the NLKP pair separation in a detector positioned 
at the same depth in the Earth as the IceCube detec- 
tor (centered at 1950 m deep); the reduced scattering in 
the lower-Z ice is neglected. The separation is shown for 
three values of the NLKP mass. The separation distri- 
butions are based on 30,000 simulated events. The sep- 
aration due to the initial opening angle is described in 
The multiple scatter enhancement to the separation 
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FIG. 1: Distribution of 300 GeV NLKP track separation 
distance in the detector. The separation is shown with and 
without multiple scattering as labeled and both cases have 
the same number of events. 



is described in the previous section. 

Multiple scattering has the largest effect on the separa- 
tion distance for NLKPs with a 300 GeV mass. The num- 
ber of events per km^ per year, at a depth equivalent to 
the IceCube detector with more than 100 m separation, 
rises from 137 to 200 events when multiple scattering is 
included, a 46% increase. For 600 and 900 GeV NLKPs, 
the corresponding increases are 30% and of 21% respec- 
tively. The background will be separated at most by 100 
meters, as can be seen in Figure 7 of As described in 
the next section, 100 m is the average detectable separa- 
tion in km^ telescopes. 

NLSPs are produced from neutrino interactions in the 
Earth. This interaction produces a charged left-handed 
slepton and a squark as described in Ref. [H H] . These 
particles decay immediately producing the two 150 GeV 
NLSPs. Here we show results for a 300 GeV squark. Mul- 
tiple scattering has a similar effect on 600 and 900 GeV 
squarks. 

Figure [3] shows the track separation between two NL- 
SPs going through the detector. The distributions show 
the separation due to the initial opening angle of the pair 
with and without multiple scattering. The separation 
distance between the two NLSP is the main observable 
to discriminate it from the dimuon background. 

The separation distributions shown in Figure [3] are 
based on 30,000 simulated events as described in Ref. 
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FIG. 2: Same as Figure [T] but for 600 and 900 GeV NLKPs. 



[2]. The darker (blue) histogram represents the separa- 
tion without multiple scattering while the lighter (pink) 
the separation with multiple scattering. 

The multiple scatter effect in this case is weaker than 
for NLKPs. As the NLSP is lighter than NLKPs, they 
are produced closer to the detector than NLKPs. The 
number of NLSP pairs with more than 100 meter sepa- 
ration will increase by 3%. However multiple scattering 



4 




track separation (m) 



FIG. 3: Track separation distance for NLSP pair events, in 
scenarios with a 300 GeV/c^ squark mass, with and without 
multiple scattering. The two curves have the same number 
of events. The multiple scattering enhancement on the sep- 
aration distance is about the same for 600 and 900 GeV/c^ 
squark masses. 



enhances the detection of vertical upcoming events as 
discussed in the next section. 



IV. DETECTION 

Next generation neutrino telescopes, hke IceCube [13] 
and kmSnet [H] will be able to detect these pairs of 
upward-going charged tracks. The NLSPs/NLKPs pro- 
duction was determined always assuming they will be 
produced below the horizon in relation to the detector. 
This avoids the large background of multiple muons pro- 
duced in cosmic-ray air showers. High px muons in air 
showers could mimic widely separated pairs. Isolated 
downgoing muons have been observed as far as 400 m 
from a shower core jl5j . 

The minimum detectable separation for charged pairs 
is an experimental question. However, a reasonable guess 
is that the minimum separation for nearly vertical pairs is 
comparable to the string spacing, 125 m in IceCube. For 
pairs that are more horizontal, the smaller separations 
may be visible, since, in the vertical direction, optical 
modules are spaced every 17 m, providing finer pixella- 
tion. Here, we will assume that the average minimum 
separation distance is 100 m. In kmSnet, the spacing is 



unknown, but this is still a reasonable figure. 

A detectable separation of 100 m would eliminate al- 
most completely the di-muon background. This can be 
seen in Figure 7 of Ref . [2\ . Multiple scattering will how- 
ever enhance NLSP/NLKP detection since it enhances 
the separation of events near the 100 m threshold for 
detecting two tracks. 

As described in the last section, with the 100 me- 
ter minimum distance, multiple scattering increases the 
number of visible 300 GeV (600 and 900 GeV) NLKPs 
by 46% (30 and 21%). NLSPs will have their number 
increased by 3%. Another feature to be accounted for is 
the detection dependence on the zenith angle. 

The pair separation also depends on the charged par- 
ticle arrival direction. Vertical upcoming particles can 
be produced farther from the detector. Therefore they 
have larger separation distances in the detector and also 
the multiple scattering effect is larger when compared to 
horizontal events. 

Multiple scattering increases the number of near ver- 
tical upcoming NLSPs pairs separated by more than 
100 meters by 24%. Here, near vertical(horizontal) is 
defined as pairs within a 5 degrees of vertical (horizon- 
tal) direction. For the near horizontal events multiple 
scattering has an insignificant effect the increase in the 
number of NLSPs separated by more then 100 meters is 
only 1.8%. Although vertical upgoing NLSPs are a tiny 
fraction of the event rate, we show the enhancement ef- 
fect to highlight the fact that multiple scattering effect 
increases with increasing zenith angle. In this way, the 
far from horizontal events will be more affected by mul- 
tiple scattering. 

For NLKP pairs multiple scattering is also less effective 
for near vertical events when compared to near horizon- 
tal events. As an example, 300 GeV near vertical pairs 
will have an 67% increase in the number of events with 
more than 100 m separation when multiple scattering is 
included while near horizontal events will have a corre- 
spondent 16% increase. 



V. CONCLUSIONS 

We have determined the multiple scattering effect on 
NLKP and NLSP particles transversing the Earth. We 
have shown that multiple scattering will enhance the sep- 
aration between NLKP pairs and between far from hori- 
zontal NLSPs. 

The pair separation in the detector is a clear signa- 
ture both for NLKPs as for NLSPs. As discussed in sec- 
tion |IV| 100 meters separation is an average lower limit 
on the experimental separation detectability. Although 
this separation by itself gets rid of most of the dimuon 
background, it will still enhance the detectability of the 
next to lightest particles. 

We find a 46% increase on the number of 300 GeV 
NLKPs separated by more than 100 meters and 30 and 
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21% for 600 and 900 GeV NLKPs when multiple scatter- 
ing is taken into account. 

The effect is less significant for NLSPs. The number 
of pairs separated by more than 100 meters has a 3% 
increase when multiple scattering is taken into account. 
Scattering can, however, enhance the number of NLSPs 
that are produced far from the horizontal direction. 

In conclusion, NLKP multiple scattering through the 
Earth is an important effect to be considered. This im- 



portance of scattering increases for pairs that are closer 
to vertical. 

This work was partially funded by the U.S. Na- 
tional Science Foundation under grant number 0653266, 
the U.S. Department of Energy under contract num- 
bers DE-AC-76SF00098 and DE-AC02-07CH11359 and 
the Brazilian National Counsel for Scientific Research 
(CNPq). 



[1] I. Albuquerque, G. Burdman and Z. Chacko, Phys. Rev. 

Lett. 92, 221802 (2004). 
[2] I. Albuquerque, G. Burdman and Z. Chacko, Phys. Rev. 

D 75, 035006 (2007). 
[3] I. Albuquerque, G. Burdman, C. Krenke, and B. Nosrat- 

pour, Phys. Rev. D 78, 015010 (2008). 
[4] B. Canadas, D. G. Cerdeno, C. Munoz and S. Panda, 

|arXiv:0812.l"067| [hep-ph] . 
[5] M. H. Reno, I. Sarcevic and S. Su, Astropart. Phys. 24, 
107 (2005) and Y. Huang, M. H. Reno, I. Sarcevic and 
J. Uscinski, Phys. Rev. D 74, 115009 (2006). 
[6] C. Amsler et ai, Phys. Lett. B667, 1 (2008). 
[7] Y. S. Tsai, Rev. Mod. Phys. 46, 815 (1974). 
[8] R. Gandhi, C. Quigg, M. H. Reno and L Sarcevic, As- 
tropart. Phys. 5, 81 (1996J 

[9] http://hyperphysics.phy-astr.gsu.edu/Hbase/tables/elabund.html' 
[10] http: / /www. mpch-mainz.mpg.de/~jesnow/Ozeanbodcn/2001 /Lecture4/Composition.html 
[11] http://en.wikipedia.Org/wiki/Earth#Chemical_compositionj 
[12] N. Olsen et ai, Geophys. J. Int. 166, 67 (2006). 
[13] S. Klein, preprint arXiv:0807.0034 
[14] C. Markou, Nucl. lustrum. & Meth. A595, 54 (2008). 
[15] S. Klein and D. Chirkin, in preprint arXiv:0711.0353 



